First— principles calculation of the intersublattice exchange interactions and Curie 
temperatures of full Heusler alloys Ni 2 MnX (X=Ga, In, Sn, Sb) 



E. §a§ioglu, L. M. Sandratskii and P. Bruno 
Max-Planck Institut fur Mikros trukturphysik, D-06120 Halle, Germany 
(Dated: February 2, 2008) 

The interatomic exchange interactions and Curie temperatures in Ni-based full Heusler alloys 
Ni2MnX with X=Ga, In, Sn and Sb are studied within the framework of the density- functional the- 
ory. The calculation of the exchange parameters is based on the frozen-magnon approach. Despite 
closeness of the experimental Curie temperatures for all four systems their magnetism appeared to 
differ strongly. This difference involves both the Mn-Mn and Mn-Ni exchange interactions. The 
Curie temperatures, To, are calculated within the mean-field approximation by solving a matrix 
equation for a multi-sublattice system. Good agreement with experiment for all four systems is 
obtained. The role of different exchange interactions in the formation of Tc of the systems is 
discussed. 
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I. INTRODUCTION 

Much efforts is currently devoted to the search for fer- 
romagnetic materials suitable for application in semicon- 
ductor spintronics devices. To allow an efficient spin- 
injection into semiconductor these materials must sat- 
isfy a number of conditions. In particular they must 
have the Curie temperature noticeably higher than the 
room temperature, be compatible with the semiconduc- 
tors used industrially, and possess a very high spin- 
polarization of the electron states at the Fermi level. 1 
Some of the Heusler compounds were found to have half- 
metallic ground stated which is characterized by a 100% 
carrier spin-polarization. 

A feature of other Heusler alloys that is useful for 
spintronic applications is very small lattice mismatch 
with widely employed semiconductors (e.g., Ni2MnIn and 
InAs)A^ Among further prominent physical properties 
of this class of materials one can mention a martensitic 
transformation in N^MnGa which takes place below the 
Curie temperature^ An interesting combination of phys- 
ical properties makes Heusler alloys the subject of in- 
tensive experimental and theoretical investigations which 
could motivate their use as multifunctional materials in 
practical applicationsi 8 ' 9 ' 10 : 11 ' 12 : 13 ' 14 

In the present work we report the theoretical study of 
the exchange interactions and Curie temperature of four 
foil Heusler compounds N12M11X with X=Ga, In, Sn and 
Sb. Experimentally, all of them are ferromagnetic and 
have similar values of the Curie temperature^ 

First important contribution to the density functional 
theory study of these systems was made in early paper 
by Kiibler et al.^ where the microscopic mechanisms 
of the magnetism of Heusler alloys was discussed on the 
basis of the calculation of the ferromagnetic and antifer- 
romagnetic alignments of the Mn moments. Recently a 
detailed study of the magnetic interactions in N^MnGa 
and N12M11AI was reported by Enkovaara et a! 16 . The 
authors used the plane spiral configurations and have 
shown that the Ni sublattice plays important role in the 



magnetic properties of the system. 

The main purpose of the present work is a detailed 
study of the exchange interactions in these systems. In 
particular we report a systematic study of the exchange 
interaction between atoms of different sublattices and 
show that pattern of exchange interactions in these sys- 
tems deviates strongly from the physical picture that can 
be expected on the basis of the experimental informa- 
tion available. Indeed common crystal structure, simi- 
lar chemical composition and close experimental values 
of the Curie temperature make the assumption natural 
that the exchange interactions in these systems are sim- 
ilar. Our study shows, however, that this assumption 
is not correct. The exchange interactions vary strongly 
depending on the X constituent. In particular, the inter- 
sublattice interactions change strongly from system to 
system. We show that different exchange interactions 
lead, in agreement with experiment, to similar values of 
the Curie temperatures. We analyze the relation between 
the properties of the exchange interactions and the Curie 
temperatures. 

The paper is organized as follows. In Sec. II we present 
the calculational approach. Section III contains the re- 
sults of the calculations and discussion. Section IV gives 
the conclusions. 



II. CALCULATIONAL APPROACH 

The calculations are carried out with the augmented 
spherical waves methodic within the atomic-sphere 
approximation^ The exchange-correlation potential is 
chosen in the generalized gradient approximation^ A 
dense Brillouin zone (BZ) sampling 30 x 30 x 30 is used. 
In all calculations the experimental values of the lattice 
parameters are used (Table QJ. 20 The radii of all atomic 
spheres are chosen equal. 

We describe the interatomic exchange interactions in 
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terms of the classical Heisenberg Hamiltonian 

= £ J RR S R*R' (i) 

R,R' 
(/iR^i/R') 

In Eq.JQ), the indices /x and v number different sublat- 
tices and R and R' are the lattice vectors specifying the 
atoms within sublattices, is the unit vector pointing 
in the direction of the magnetic moment at site (/x,R). 
The systems considered contain three 3d atoms in the 
unit cell with positions (2,7,7) for the Mn atom and 
(0, 0, 0) and (5,5,5) for two Ni atoms. 

We employ the frozen-magnon approach2ii22i2ii to cal- 
culate interatomic Heisenberg exchange parameters. The 
calculations involve few steps. In the first step, the ex- 
change parameters between the atoms of a given sublat- 
tice n are computed. The calculation is based on the 
evaluation of the energy of the frozen-magnon configura- 
tions defined by the following atomic polar and azimuthal 
angles 

0R = 0> R -q-R+0 M - (2) 

In the calculation discussed in this paper the constant 
phase </) M is always chosen equal to zero. The magnetic 
moments of all other sublattices are kept parallel to the 
z axis. Within the Heisenberg model (JTJ the energy of 
such configuration takes the form 

£?w (0, q) = (9) + sin 2 6J^ (q) (3) 

where Eq m does not depend on q and the Fourier trans- 
form J Miy (q) is defined by 

J^(q)=£j^exp(iq-R). (4) 

R 

In the case of v = fi the sum in Eq. J2J does not in- 
clude R = 0. Calculating E^^(9,q) for a regular q- 
mesh in the Brillouin zone of the crystal and performing 
back Fourier transformation one gets exchange parame- 
ters Jq^ for sublattice \i. 

To determine the exchange interactions between the 
atoms of two different sublattices fi and v the frozen- 
magnon configurations (Eq. [2J are formed for both sub- 
lattices. The Heisenberg energy of such configurations 
takes the form 

E^(9,q) = E^(0)+sm 2 0[J^{q) + J""(q)] 

+2sin 2 <9ReJ^(q) (5) 

where Eq U (9) is a q-independent part. Performing calcu- 
lation of [E^(6,q) - E^(6,0)} and subtracting single- 
sublattice contributions known from the previous step 
one finds [ ReJ^(q) - ReJ^(O)]. The back Fourier 
transformation of this expression gives for R ^ the 
following combinations of the interatomic exchange pa- 
rameters: = 5 (-'or, + ^o(-R))' I n S enerai , one needs 



to perform similar calculations for different phases </> M 
in Eq. @ to get another linear combination of the pa- 
rameters to be able to separate them. For the systems 
considered these calculations can be avoided by taking 
into account the symmetry of the lattice (see below). 

Quantities does not contain information about the 
interaction of the atoms within the first unit cell corre- 
sponding to R = 0. These exchange parameters can be 
found in the calculations for magnetic configurations pe- 
riodic with the periodicity of the lattice. The atoms in 
the unit cell are separated into two groups. Within each 
group the moments are parallel. The moments from the 
different groups form an angle 8. The energies of such 
magnetic configurations are expressed through the sums 
Jq V = J2r ^or- Since the sums X^r/o ^or are known 
from the preceding step the parameters with R = be- 
come accessible. The symmetry relation Jqq — Jq R for 
R = (0, 5, 5) allows to split the sums to individual 
parameters. 

The Curie temperature is estimated within the mean- 
field approximation for a multi sublattice material by 
solving the system of coupled equations 24 

where (s u ) is the average z component of sj^. Eq. © can 
be represented in the form of eigenvalue matrix problem 

(© - TI)S = (7) 

where 0^„ — j^Jq" , I is a unit matrix and S is the 
vector of (s"). The largest eigenvalue of matrix O gives 
the value of Curie temperature^ 

III. RESULTS AND DISCUSSION 

In Table [I] we present calculated magnetic moments. 
For comparison, the available experimental values of the 
moments and the results of previous calculations are pre- 
sented. Relative variation of the Mn moment is small. 
On the other hand, the moment of Ni and X atoms show 
strong relative variation and are in N^MnSb about two 
times smaller than in Ni 2 MnGa or Ni 2 MnIn. The values 
of the magnetic moments are in good agrement with the 
results of previous calculations. 

The calculated Heisenberg exchange parameters are 
presented in figure ^ As mentioned in the introduction, 
the assumption that the closeness of the experimental 
Curie temperatures is the consequence of the similarity 
of the exchange interactions is not confirmed by the cal- 
culations. We obtain strong dependence of the exchange 
interactions on the type of the X atom. For X=Ga and 
X=In that belong to the same column of the Mendellev's 
table (see inset in Fig. [2J) we obtain similar pattern of 
Heisenberg exchange parameters. On the other hand, 
for X atoms belonging to different columns the changes 
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TABLE I: Experimental lattice parameters and magnetic mo- 
ments (in fi B ) of Ni 2 MnX (X=Ga, In, Sn, Sb). 





a(a.u.) 


Mn 


Ni 


X 


Cell 


Ni 2 MnGa 


11.058" 


3.570 


0.294 


-0.068 


4.090 






(3.43 6 ) 


(0.36 6 ) 


(-0.04 6 ) 


(4.11 6 ) 












(4.17 c ) 


Ni 2 MnIn 


11.468° 


3.719 


0.277 


-0.066 


4.208 


Ni 2 MnSn 


11.419 a 


3.724 


0.206 


-0.057 


4.080 






(3.53 d ) 


(0.24 d ) 


(-0.03 d ) 


(4.08 d ) 


Ni 2 MnSb 


11.345" 


3.696 


0.143 


-0.033 


3.950 



a Rcf 20 

& Ref|2a 
c Ref|U (Exp.) 
d Ref|23 
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FIG. 1: The parameters of the Mn-Mn (intra-sublattice) ex- 
change interactions and Mn-Ni (inter-sublattice) exchange 
interactions in Ni 2 MnX (X=Ga, In, Sn, Sb). The distances 
are given in the units of the lattice constant. The significance 
of the oscillations of the exchange parameters is verified by 
varying the q mesh in the frozen-magnon calculations. 

in the exchange interactions are very strong (Fig. QJ. 
These changes concern both the Mn-Mn intra-sublattice 
interactions and the Ni-Mn inter-sublattice interaction. 

Considering the Mn-Mn interactions we notice that in 
Ni2MnGa and N^Mnln the interaction with the coordi- 
nation spheres from the first to the forth is positive. The 
interaction with the first coordination sphere is weaker 
than with the following ones. The interaction with the 
fifth sphere is very small. The interaction with the 6th 
sphere is negative. The interaction with further coordi- 
nation spheres is very weak. 

In Ni2MnSn the interaction with the first sphere 
strongly increases compared with N^MnGa and 
Ni2MnIn. On the other hand, the interaction with the 
third sphere becomes small. The interaction with the 
forth sphere is strongly negative. The interaction with 
further neighbors are weak. 

The trend observed in transition from Ni2MnGa and 



FIG. 2: Frozen-magnon energies as a function of the wave 
vector q (in units of 2n/a) in Ni 2 MnX for Mn-Mn (left) and 
Mn-Ni interactions (right). 



Ni2MnIn to Ni2MnSn becomes even stronger in the case 
of X=Sb. Here the interaction with the first neighbor 
increases further and is the only strong exchange inter- 
action between the Mn atoms. 

The inter-sublattice Mn-Ni interaction behaves very 
differently. A sizable interaction takes place only between 
nearest neighbors. This interaction is very strong in 
Ni2MnGa and Ni2MnIn and quickly decreases for X=Sn 
and, especially, X=Sb. 

To reveal the physical origin of the strong difference in 
the exchange parameters of these systems we plot in Fig. 
|2]the frozen-magnon energies as a function of wave vector 
q for one direction in the Brillouin zone^ We remind 
the reader that the inter-atomic exchange parameters are 
Fourier transforms of the E(q) functions and therefore 
reflect their form. Indeed, E(q) curves for Ni2MnGa and 
Ni2MnIn are close to each other that leads to a similar 
set of interatomic exchange parameters (Fig. These 
curves deviate strongly from a simple cosinusoid having 
a maximum at q about 0.6 and a rather weak variation 



NLMnGa 



Ni 2 MnIn 




Ni,MnSb 



\ 






— Maj. 
Min. 


\ > 

\ 

\ I 


1 
\ 




l 


\ 
I 

\ 1 
v\ 1 

\\v 

\ \ 





FIG. 3: Band structure of Ni2MnX along the high symmetry 
line (r-X). 



at q > 0.6i2S The complexity of E(q) means that several 
Fourier components are needed to describe the features 
of the function. This is reflected in the the Heisenberg's 
parameters of NLMnGa and N^Mnln. 

On the other hand, the E(q) curve of Ni2MnSb is well 
described by one cosinusoid (Fig. |2J) that results in a sin- 
gle large Mn-Mn exchange parameter (Fig.^J. The E(q) 
of N^MnSn assumes an intermediate position from the 
viewpoint of the complexity of the function. This prop- 
erty is also reflected in the exchange parameters (Fig-QJ. 

Note that the character of the q-dependence of the 
energy illustrated by Fig. [2 is a consequence of the prop- 
erties of the electronic structure of the compounds. In- 
deed, in Fig. |3| we see that the electronic structures of 
NLMnGa and N^Mnln are similar. Transition along the 
row In-Sn-Sb leads to increasing difference in the elec- 
tron spectrum. This increasing difference can be traced 
back to the change in the number of valence electrons: a 
Sb atoms has two more valence electrons than In and one 
more electron than Sn. As the result an important differ- 
ence in the electron structure of the system is a relative 
shift of the Fermi level to a higher energy position in the 
sequence In-Sn-Sb. This shift is clearly seen in the DOS 
presented in Fig.Q] The positions of the same features of 
the DOS in different systems are well described by linear 
functions with negative angle coefficients. For the Mn 
peaks all three lines are almost parallel. This means that 
the change in the Mn DOS from system to system can 
be treated as a rigid shift with respect to the Fermi level. 
In the case of the Ni-DOS the situation is more compli- 
cated since, besides the variation of the electron number, 
an additional influence on the peak positions is exerted 
by the variation of the Ni atomic moment (see Table QJ. 

The E(q) curves determining the Mn-Ni interactions 
are presented in Fig. [5J The form of the curves is in all 
cases close to a cosinusoid. Therefore, only one exchange 
parameter has sizable value. The strength of the inter- 
action is in correlation with the value of the magnetic 




FIG. 4: Spin-projected partial density of states of Ni2MnX. 
The separate graph in the left-hand bottom part of the figure 
shows the variation of the positions of a number of the DOS 
peaks for the In-Sn-Sb series of compounds. The lines in this 
graph are guides for the eye. The energies are measured with 
respect to the Fermi level of the corresponding system. The 
small open circles on the DOS curves mark the positions of 
the peaks. Arrows show spin projections. The analysis of the 
states at the Fermi level shows that the main contributions 
comes from the 3d states of Mn and Ni. 



TABLE II: Mean-field estimation of the Curie temperatures 
for Ni2MnX (X=Ga, In, Sn, Sb). The experimantal Curie 
temperature values are taken from Refl20l 



T, 



Mn-Mn[MFA] 



(K) 



-,[MFA] 



(K) 



,[Exp] 



Ni 2 MnGa 
Ni 2 MnIn 
Ni 2 MnSn 
Ni 2 MnSb 



302 
244 
323 
343 



389 
343 
358 
352 



380 
315 
360 
365 



moment of the Ni atoms. 

The interatomic exchange parameters are used to eval- 
uate the Curie temperature. If only the Mn-Mn interac- 
tions are taken into account we obtain values shown in 
Table [H] Despite very strong difference in the Mn-Mn 
exchange parameters in these systems the difference in 
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the corresponding Curie temperatures is not very large. 
The explanation for this result is the property that in 
MFA to a one-sublattice ferromagnet the value of the 
Curie temperature is determined by the sum of the inter- 
atomic exchange interactions Jo = X^r^o ^or- Jo gives 
the average value of E(q) and is less sensitive to the de- 
tailed form of the E(q) function. 

The comparison of the Curie temperatures calculated 
with the use of the Mn-Mn exchange parameters only 
with experimental Curie temperatures shows that the 
agreement with experiment is not in general good. In 
the case of Ni 2 MnGa the error is about 30%. 

Account for inter-sublattice interactions improves the 
agreement with experimental Tc values considerably (Ta- 
ble [HJ . This shows that the Ni moment provides a mag- 
netic degree of freedom which plays important role in the 
thermodynamics of the system. 



Heusler alloys Ni 2 MnX (X= Ga, In, Sn, Sb) within the 
parameter-free density functional theory. Our calcula- 
tions show that despite similarity of the Curie tempera- 
tures of these systems there is strong difference in the un- 
derlying magnetic interactions. This difference involves 
both the Mn-Mn and Mn-Ni exchange interactions which 
depend strongly on the X constituent. The Curie temper- 
atures are calculated within the mean-field approxima- 
tion to the classical Heisenberg Hamiltonian by solving 
a matrix equation for a multi-sublattice system. Good 
agreement with experiment for all four systems is ob- 
tained. 



IV. CONCLUSION 
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